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ABSTRACT 

The statistics of wide-separation (6" < 9 < 15") gravitational lenses constrain the amount of mass 
in the cores of dark matter halos on group and cluster mass scales. For a family of halo models with 
a central cusp p cx r~ a (1.0 < a < 1.9), the lack of wide-separation lenses in the large Cosmic Lens 
All-Sky Survey yields an upper limit on the fraction of the halo mass that is contained within of 
the virial radius, / CO re < 0.023 (95% confidence level, LCDM). This limit offers an important test of the 
cold dark matter paradigm. While the halo profiles derived from numerical simulations appear to be 
consistent with this upper limit, larger surveys currently underway such as the 2dF and SDSS should 
detect wide-separation lenses and thus provide a measurement of the core mass fraction in massive dark 
matter halos. 

Subject headings: cosmology: dark matter - gravitational lensing 



1. INTRODUCTION 

Cuspy dark matter halos are a robust prediction of the 
' cold dark matter (CDM) paradigm. Numerical simula- 
tions have shown that there is a "universal" density profile 
for dark matter halos with the form p a r" 1 (r s + r) -2 , 
where r s is a scale length (Navarro, Frenk & White 1997, 
hereafter NFW) . This result does not depend on particular 
; cosmogonies or on specific initial conditions or formation 
histories (Huss, Jain & Steinmetz 1999a, b). There is some 
disagreement about the very inner regions of halos; recent 
simulations have suggested that the central cusp may be 
| steeper than r -1 (e.g., Moore et al. 1998, 1999) and may 
■ depend on the halo mass (Jing & Suto 2000). 

It is important to compare the predicted halo pro- 
files with observational data to test the CDM paradigm. 
Galaxy rotation curves can be used to probe dark matter 
halos (e.g., Flores & Primack 1994; McGaugh & de Blok 

1998) , but the implications are not clear because beam 
smearing can make HI rotation curves appear shallower 
than they really are (van den Bosch & Swaters 2000). 

i Gravitational lensing offers an additional test of halo pro- 
files. More than 60 multiply-imaged quasars and radio 
sources are known (see Falco et al. 1999), and with image 
separations of a few arcseconds they primarily probe halos 
on galaxy mass scales. These lenses can be used to con- 
strain CDM halos (e.g., Rusin & Ma 2000), but the test 
is complicated by the necessity of including baryons and 
their effects on dark matter halos. At least 24 giant lensed 
arcs are also known (see Williams, Navarro & Bartelmann 

1999) , and with image separations of several tens of arc- 
seconds or larger they probe rich clusters. 

Lenses with intermediate image separations, 9 ~ 10", 
probe halos at intermediate masses. There is just one con- 
firmed lens (Q 0957+561, Walsh, Carswell & Weymann 
1979) and one good candidate (RX J0921+4529, Muhoz et 
al. 2000) with 9 > 6", and each is produced by a galaxy in a 
cluster. No systematic lens survey has found any lenses in 
this regime, and this result has been used to place limits on 
cosmological parameters (e.g., Cen et al. 1994; Kochanek 
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1995; Wambsganss et al. 1995; Maoz et al. 1997; Mortlock 
& Webster 2000). Instead, in this Letter we propose to 
use the statistics of wide-separation lenses as a new test of 
the structure of dark matter halos (see Flores & Primack 
1996). We consider a range of halo profiles suggested by 
recent CDM simulations, and compare our model predic- 
tions with the latest and largest lens survey, the Cosmic 
Lens All-Sky Survey. 

2. CUSPY HALOS 

To obtain a family of halos with a range of central cusps, 
we adopt a generalization of the NFW profile (see Jing & 
Suto 2000; Wyithe, Turner & Spergel 2000), 
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where r s is a scale length and p s is a characteristic density. 



The density scales as p cx r~ a for small r and p cx r~ 
large r, so it reduces to the NFW profile when a = 1. 
characteristic density p s is given by 
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where 2-Fi(a, b, c, z) is a hypergeometric function, and f2oo 
is the radius within which the mean density is 200 times 
the critical density of the universe p C rit- Moore et al. (1998, 
1999) advocate a profile of the form p cx x~ 15 (l +X 1 ' 5 )" 1 
where x — r/r s . Although this differs slightly from a cuspy 
profile with a — 1.5, the two profiles are similar at small 
radii and differ by <;15% in the enclosed inner mass, which 
we will argue below is the most important quantity in de- 
termining the number of lenses. We do not consider the 
Moore profile further because there is no simple family of 
models connecting it to the NFW model. It is often con- 
venient to replace the scale radius of a cuspy halo with 
a "concentration" parameter. For NFW (a = 1) halos, 
Navarro et al. (1997) define Cnfw = ^200 As- For the gen- 
eralization to a 7^ 1, we define C_2 = ?'20o/ r -2 in terms of 
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the radius r_2 at which the logarithmic slope of the density 
profile is —2, which is equivalent to Cnfw f° r a — I- 3 

A cuspy halo is fully specified by its mass, redshift, cusp 
slope, and scale radius or concentration, but these param- 
eters are not all independent. For NFW halos, the con- 
centration is correlated with the mass and redshift in a 
way that appears to contain information about the forma- 
tion epoch of the halo (e.g., Navarro et al. 1997). How- 
ever, there is a scatter in concentrations at fixed mass and 
redshift that may reflect differences in formation histo- 
ries (e.g., Bullock et al. 1999; Klypin et al. 2000). The 
scatter is important for our calculations because lens- 
ing is very sensitive to concentration; more concentrated 
halos are much more efficient lenses. We use results 
from the simulations by Bullock et al. (1999), namely a 
scatter that is consistent with a log-normal distribution 
with standard deviation A(logC_2) — 0.18, and a me- 
dian concentration that scales with mass and redshift as 
med(C_ 2 ) oc M-^l + z)" 1 (Bullock et al. 1999; Navarro 
& Steinmetz 2000). These trends were derived from fits 
of NFW (a = 1) profiles to simulated halos, but because 
their origin is physical we assume that they apply to other 
values of a as well. 

3. LENS THEORY 

To compute the total number of lenses expected to 
be produced by a population of halos it is sufficient to 
use spherical halos. Departures from spherical symmetry 
mainly affect the relative numbers of two-image and four- 
image lenses, which we do not differentiate. Wyithe et al. 
(2000) and Li & Ostriker (2000) discuss the lensing prop- 
erties of cuspy halos of the form given in eq. ([j]) , but their 
formalism can be simplified so the deflection is written as 
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where M{r) is the mass inside radius r and S cr = 
(c 2 D s ) / (4itGDiDi s ) is the critical surface density for lens- 
ing (e.g., Schneider, Ehlers & Falco 1992). Here Di and D s 
are angular diameter distances to the lens and source, re- 
spectively, and Di s is the angular diameter distance from 
the lens to the source. The positions of lensed images 
are found by solving the lens equation u = R — (f>^ R (R), 
where u is the impact parameter of the source relative 
to the lens. The magnification of an image at R is then 
At = [1 - (f>, R /R]- l [l - ^.flfl] -1 , where <f) :RR = d(<f) :R )/dR. 

The number of lenses expected to be found in a survey 
depends on the optical depth for lensing as well as on the 
"magnification bias," which accounts for lenses that are in- 
trinsically fainter than the flux limit of the survey but are 
brought into the sample by the magnification (e.g., Turner, 
Ostriker & Gott 1984). The number of lenses with a total 
flux greater than 5* expected to be found in a survey with 
selection functions described by T is then 
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3 For a^l, C—2 differs from the parameter C s = T20o/t s used by 
q/1 because the radius r_2 has physical significance, while r s does 
halos, C s = (2 - a)C- 2 - 



where z s is the source redshift, dV is the comoving volume 
element, and dn/dM is the mass function of halos. Also, u 
is the angular position of the source relative to the lens cen- 
ter, and fj, is the total magnification of that source; the u 
integral extends over the range of impact parameters that 
produce multiple images. Finally, [dN src (> S)/dz s ] dz s is 
the number of sources brighter than flux S that lie in the 
redshift range z s to z s + dz s . The factor T{u) indicates 
whether a lens associated with a source at u would be 
detected given the selection functions. 

We follow Narayan & White (1988), Kochanek (1995), 
and Porciani & Madau (2000) and compute lens statis- 
tics using a mass function of dark matter halos given 
by Press-Schechter theory combined with the spherical 
collapse model. Numerical simulations suggest that this 
mass function overestimates the number of halos below 
~ 10 14 hr 1 Mq and underestimates the number of halos 
above this mass (e.g., Jenkins et al. 2000). Because we 
are interested in the high-mass end of the mass func- 
tion, our results should slightly underestimate the num- 
ber of lenses, and lead to conclusions that are conserva- 
tive. We compute the CDM power spectrum using the 
fitting formula given by Eisenstein & Hu (1999). We 
present results for two different flat cosmologies: SCDM, 
with matter density Qm = 1 and Hubble constant Hq = 
50 km s" 1 Mpc' 1 ; and LCDM, with Q M = 0.3, 0\ = 0.7, 
and Hg = 65 km s _1 Mpc -1 . From the abundance of clus- 
ters we set the variance of mass fluctuations on8/i _1 Mpc 
scales to be a 8 = 0.52 for SCDM and ct 8 = 0.93 for LCDM 
(Eke, Cole & Frenk 1996). 

Li & Ostriker (2000) have recently performed a similar 
calculation. However, they use a different power spectrum, 
omit magnification bias, and do not consider a scatter in 
concentrations, so their quantitative results are somewhat 
different from ours. 



4. THE CLASS LENS SURVEY 

To test the models we must compare predicted lens 
statistics with a well-defined observational sample. The 
Cosmic Lens All-Sky Survey (CLASS; e.g., Browne & My- 
ers 2000) is the largest statistically homogeneous search 
for gravitational lenses. The sample comprises 10,499 
flat-spectrum radio sources whose flux distribution can be 
described as a power law dN SIC /dS oc S* -21 (see Rusin 
& Tegmark 2000). The survey includes 18 gravitational 
lenses with image separations 9 < 3". An explicit search 
for lenses with image separations 6" < 6 < 15" has found 
no lenses (Phillips et al. 2000). 

The redshift distribution of the full CLASS sample is 
not known. Marlow et al. (2000) report redshifts for a 
small subsample of 27 sources. They find a mean redshift 
of (z) = 1.27, which is comparable to that found in other 
radio surveys at comparable fluxes. We present results as- 
suming that the full sample has the same source redshift 
distribution as the spectroscopic subsample of Marlow et 
al. (2000), but our results are not substantially different 
if we simply place all sources at the mean redshift of the 
spectroscopic subsample, so the source redshift distribu- 
tion is not a significant source of uncertainty. 

Wyithe et al. (2000). We believe that C—2 is the better generalization to 
not; in any case, the alternate parameters have a simple relation in cuspy 
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5. RESULTS 

We focus on wide-separation lenses, 6" < 9 < 15", which 
are not present in the CLASS sample and correspond to 
halos on the mass scale of groups and clusters. The exact 
bounds on the image separation have no intrinsic signifi- 
cance but are chosen to correspond to the analysis of the 
CLASS survey data (see Phillips et al. 2000). We com- 
pute the number of lenses for particular values of the cusp 
slope and median concentration; we consider a wide range 
of models with 1.0 < a < 1.9 and 2 < med(C_ 2 ) < H. 
Models with med(C_2) > 11 arc ruled out at 99% confi- 
dence or better for all values of a that we consider. 
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cusp slope a 

Fig. 1. — The number of wide-separation lenses expected in the 
CLASS survey as a function of the cusp slope and concentration 
parameter, for the LCDM cosmology. The dotted curves show con- 
tours drawn at 1, 2, 3, and 4 lenses, and the solid curves show the 
1<t, 90%, and 95% confidence upper limits. 

Figure 1 shows the expected number of wide-separation 
lenses as a function of a and C_2- The results are not easy 
to interpret for two reasons. First, there is some ambiguity 
in the meaning of the concentration parameter when ha- 
los are allowed to have arbitrary cusps (see §5 of Wyithc 
et al. 2000 for a detailed discussion). Second, lensing is 
not very sensitive to the cusp slope and the concentration 
separately, but depends mainly on the amount of mass in 
the central regions of halos. It is more instructive then to 
plot the number of lenses versus the "core mass fraction" 
/core) or the fraction of the total mass of the halo that 
is contained within some small fiducial radius, as in Fig- 
ure 2. The best fiducial radius is the one that minimizes 
the scatter in the relation between the core mass fraction 
and the number of lenses. 4 It is 0.045 r2oo for SCDM and 
0.032 r2oo for LCDM, which is approximately equal to the 
mean Einstein radius of lenses with 6" < 9 < 15". Because 
of systematic trends and the scatter in concentrations, the 
core mass fraction is measured for a median halo of mass 
10 14 h~ x Mq at redshift z = 0, which is the mass scale 
relevant for wide-separation lenses. 
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Fig. 2. — The number of wide-separation lenses plotted versus 
the fraction of halo mass that is contained within a small fidu- 
cial radius (indicated in the label). Each point indicates a model 
with specific values of a and med(C— 2), with squares (triangles) for 
LCDM (SCDM). The filled circles indicate NFW (a = 1) models 
with med(C_2) = 6, 7, and 8. Lines connect models with a = 1 
and 1.5 (in the main relations), and a = 1.8 and 1.9 (offset above 
the main relations). The horizontal lines indicate the la, 90%, and 
95% confidence upper limits on the number of lenses. 

The number of predicted lenses is strongly correlated 
with the core mass fraction, and this relation explains 
most of the dependence of the number of lenses on the 
cusp slope and the concentration. This result is surprising 
in its simplicity, because lensing actually depends on the 
projected mass distribution and because it is not obvious 
that lensing should measure the fraction of mass that is 
contained within a fiducial radius expressed as a fraction 
of the halo's radius. Indeed, for these reasons the relation 
between lensing and the core mass fraction is not perfect; 
models with steep cusps (a — 1.8 or 1.9) fall somewhat 
off the main relation, and softened isothermal models (not 
shown) follow a similar but slightly different relation. Nev- 
ertheless, the main relation between the number of lenses 
and the core mass fraction is quite tight for a range of cusp 
slopes and concentrations that is wider than the range in- 
dicated by numerical simulations. Our first result, then, is 
that the statistics of wide-separation lenses offer a surpris- 
ingly good method for constraining the fraction of mass in 
group and cluster halos contained within a fiducial radius 
that is a few percent of r 2 oo- 

The predicted number of lenses admits the following 
quantitative interpretation. The probability of finding 
k lenses when N are expected is well approximated by 
the Poisson distribution as P(k\N) = N k e~ N /k\. Given 
that CLASS contains no wide-separation lenses (k = 0), 
the upper limits on the predicted number of lenses are 
N < (1.15, 2.30, 3.00) at the (la, 90%, 95%) confidence lev- 
els. Figures 1 and 2 then yield upper limits on the core 
mass fraction, cusp slope, and concentration, which are 
summarized in Table 1. For NFW (a = 1) halos, the 
median concentration of a 10 14 hr 1 M Q halo must be less 
than about 7 at la and less than 9 at 90%. The simu- 



4 Specifically, we minimize the scatter in the log-log relation, including all models with 1.0 < a < 1.7 and 2 < mcd(C-2) < 11 that predict 
at least 0.1 lenses. 
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lations by Navarro et al. (1997) and Bullock et al. (1999) 
yield concentrations of slightly less than 6, so they are con- 
sistent with the lens data (see Figure 2). Halos that are 
significantly more concentrated would disagree with the 
lens data. We could compute the limits on the concentra- 
tion parameter for other cusp slopes (and Table 1 includes 
the limits for Moore- like a = 1.5 halos), but we believe it 
is more instructive to focus on the core mass fraction. 

Finally, we note that modeling all halos as singular 
isothermal spheres yields a prediction of 9 (11) wide- 
separation lenses for SCDM (LCDM). Thus, the statistics 
of wide-separation lenses strongly exclude the hypothesis 
that the halos of groups and clusters are singular isother- 
mal spheres. 



Upper Limits From Wide-Separation Lenses 



Cosmology 


Quantity 


la 


90% 


95% 


SCDM 


/core 


0.031 


0.037 


0.040 




C-2(l) 


7.6 


9.0 


9.7 




C_ 2 (1.5) 


5.6 


7.3 


8.1 


LCDM 


./core 


0.017 


0.021 


0.023 




C-a(l) 


7.3 


9.0 


9.8 




C_ 2 (1.5) 


4.7 


6.6 


7.5 



NOTE. — /core is the fraction of a median halo's mass that is con- 
tained within 0.045 r 20 o for SCDM, or 0.032 r 20 o for LCDM. C_ 2 (a) 
is the median concentration of halos with cusp slope a. 

6. CONCLUSIONS 

We have studied the statistics of wide-separation (6" < 
8 < 15") lenses produced by dark matter halos with gen- 
eral cusps of the form p oc r~ a for 1.0 < a < 1.9. In 
these models the number of expected lenses is determined 
almost entirely by the fraction of the halo mass that is 
contained within a fiducial radius that is ~4% of the virial 
radius. Combining our results with the lack of observed 
wide-separation lenses in the CLASS lens survey yields an 
upper limit on how concentrated halos can be. The con- 
centration parameter for the standard NFW profile must 
be C <; 7 (ler); see Table 1 for more details. Massive halos 
found in numerical simulations of CDM cosmogonies have 
C ~ 6 (e.g., Navarro et al. 1997; Bullock et al. 1999) 
and thus are consistent with the lensing constraints, but 
halos cannot be much more concentrated than that. Ha- 
los with cusps steeper than r -1 are also allowed provided 
that they have somewhat lower concentrations. Although 



we quote limits on the concentration, we believe that it is 
more attractive to state our results in terms of the core 
mass fraction, e.g., / corc < 0.023 at 95% confidence for 
LCDM. This quantity should be a more robust prediction 
of numerical simulations than the slope of the central den- 
sity and should allow a simple comparison of our models 
with simulations. 

Two points are important for interpreting out results. 
First, the scatter in halo properties requires that we quote 
median properties. Any given halo may be more concen- 
trated than the median and hence violate the mass limits 
that we have derived for median halos. In other words, 
the lensing limits must be applied statistically, not to in- 
dividual halos. Second, we have considered only dark mat- 
ter models, neglecting any baryons that might occupy the 
cores of massive dark matter halos as cD galaxies in clus- 
ters or central ellipticals in groups. This approach should 
lead us to conservative conclusions, because any central 
baryonic component would only increase the number of ex- 
pected lenses and thus aggravate the discrepancy between 
the data and the CDM models. 

While our results indicate that the lack of wide- 
separation lenses in the CLASS survey is perhaps not sur- 
prising, they also suggest that we should expect to discover 
them soon. In our models, NFW halos with a median con- 
centration C = 6 predict 0.6 lenses with 6" < 8 < 15" 
out of the 10,499 CLASS targets. Several new surveys are 
large enough that we should expect a few group and cluster 
lenses; they include 2dF and SDSS, which should measure 
the redshifts of more than 25,000 and 10 5 quasars, respec- 
tively. If the new surveys continue to lack wide-separation 
lenses, the lensing constraints on the cores of dark mat- 
ter halos will become strong enough to question the CDM 
paradigm. If, on the other hand, the new surveys do detect 
wide-separation lenses, the lenses will allow empirical mea- 
surements of the properties of the cores of massive dark 
matter halos. 

We thank M. Steinmetz and V. Eke for helpful dis- 
cussions about dark matter halos, and P. Phillips and I. 
Browne for information about the CLASS survey. Support 
for this work was provided by NASA through ATP grant 
NAG5-4236 and grant AR-06337.10-94A from the Space 
Telescope Science Institute (P.M.). 
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